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Abstract 

Point-form relativistic quantum mechanics is applied to elastic electron- 
deuteron scattering. The deuteron is modeled using relativistic interactions 
that are scattering-equivalent to the nonrelativistic Argonne vis and Reid '93 
interactions. A point-form spectator approximation (PFSA) is introduced to 
define a conserved covariant current in terms of single-nucleon form factors. 
The PFSA is shown to provide an accurate description of data up to mo- 
mentum transfers of 0.5 GeV 2 , but falls below the data at higher momentum 
transfers. Results are sensitive to the nucleon form factor parameterization 
chosen, particularly to the neutron electric form factor. 
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I. INTRODUCTION 



Electron scattering is considered to be an ideal tool to study the electromagnetic struc- 
ture of hadronic systems. Relativity cannot be ignored for momentum transfers that provide 
information about the structure of the hadrons at the scale of a few tenths of a fermi. In 
order to understand hadronic systems at this scale, consistent relativistic models of both 
the hadronic dynamics and the hadronic electromagnetic current operator are required. If 
the dynamics and the current operator satisfy cluster properties |T]J1, then the informa- 
tion learned about the structure of the simplest two- and three-body systems provides the 
essential components of models needed to treat complex targets. 

Elastic electron-deuteron scattering is the simplest reaction that must be accurately mod- 
eled in order to constrain the dynamical generators and current operators that are needed 
to model complex systems. Because the deuteron is an isoscalar target, it might be expected 
that it can be accurately described by a pure impulse approximation. Unfortunately, pure 
impulse approximations are not consistent with current covariance or current conservation. 
An important goal is to find a physically motivated extension of the impulse approxima- 
tion that is consistent with current conservation and covariance and is also qualitatively 
consistent with experiment. 

In the one-photon-exchange approximation the experimental observables can be ex- 
pressed in terms of matrix elements of the hadronic electromagnetic current operator be- 
tween the initial and final eigenstates of the hadronic Hamiltonian. The general form of 
these matrix elements is 

(p'f^\Ho)\pj H }, (i) 

where \p' f fj/j) and \p j fij) are eigenstates of the four-momentum, spin, and three- 
component of the spin, in reference frames related by a boost with momentum transfer 
Q = p' — p. 7^(0) is the hadronic current density at x = 0. 

The dynamical constraints on the current are current conservation, 
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[P M) >(0)] = 0, (2) 

and current covariance, 

U{A, a)J^{x)U\A, a) = (Ar l y u J v (kx + a). (3) 

Here U (A, a) is the unitary representation of the inhomogeneous Lorentz group, whose ex- 
istence is required by relativistic invariance ||. P M is the four- momentum operator, with 
U(I,a) = e ipa . 

In applications there are essentially two approaches used to compute the hadronic current 
matrix elements. These are the covariant and direct interaction |p|-|T7[| approaches. 

Each approach has its own advantages and disadvantages. 

Covariant approaches assume that the underlying theory is a local quantum field theory. 
For the case of the deuteron the input is a covariant current vertex of the form 

(0|T(*(x 1 )*(a; 2 )J' t (a;)*( 3 / 1 )*( 3 / 2 ))|0), (4) 

which is the vacuum expectation value of a time-ordered product of nucleon and current 



fields. Assuming the existence of an underlying quantum theory, Mandelstam [18[ showed 
how to extract the desired current matrix elements from the vertex. The Fourier transform 
of the vertex has pole terms on the deuteron mass shell. The residue includes a pair of 
Bethe-Salpeter amplitudes and a current matrix element. The Bethe-Salpeter normalization 
condition fl9| can be used with the solution of the homogeneous Bethe-Salpeter equation 
to remove the amplitudes from the residue. What remains is the desired current matrix 
element. 

Quasipotential equations j|-|J are based on these same concepts, but they introduce 
constraints designed to preserve the physical singularities. Current matrix elements are 
extracted using the constrained amplitudes and vertex functions. 

Covariant methods are appealing because of their formal connection to a quantum field 
theory; however in most applications it is necessary to model the vertex and Bethe-Salpeter 
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kernel, and to replace the two-point Green's function with the free two-point function. Quasi- 
potential methods lead to simpler calculations J7| than the full Bethe-Salpeter approach || , 
but the reductions complicate cluster properties. 

Direct interaction approaches attempt to construct consistent models of U (A, a) and 
J M (x) directly. The transformation properties require that both U(A, a) and J M (x) have an 
interaction dependence. Dirac [2(J addressed the problem of constructing U (A, a) by includ- 



ing interactions in some of the infinitesimal generators of U(A, a). He introduced the notion 
of forms of dynamics which minimize the number of interaction-dependent generators. The 
three main forms are: the instant form, where the interactions are in the Hamiltonian and 
Lorentz boost generators; the point form, where the interactions are in the four-momentum; 
and the front form, where interactions appear in the operators that generate transforma- 
tions transverse to a fixed light front (a three-dimensional hyperplane tangent to the light 
cone.) But while Dirac identified the different possibilities for putting interactions in se- 
lected generators, he did not show how to actually construct the Poincare generators with 
interactions. 

The first exact construction of Poincare generators with interactions was due to 
Bakamjian and Thomas [BTJ using Dirac's instant form. There are Bakamjian-Thomas-like 



constructions in each of the forms of dynamics ||22|| , and they are scattering equivalent. 



While explicit dynamical models of current operators are difficult to construct, consistent 
current matrix elements can be obtained by prescriptions that evaluate selected independent 
matrix elements using single nucleon currents. The remaining current matrix elements 
can then be determined by using covariance, current conservation and discrete symmetries. 
These generate the needed dynamical contributions to the current matrix elements. 

Direct interaction approaches provide an exact treatment of the symmetries associated 
with special relativity, but are not directly related to an underlying field theory. A number of 
direct interaction applications to elastic electron-deuteron scattering exist in the literature. 
To date most applications have used Dirac's instant- [[17] or front-forms [2B|, r5| , P^ -Pm of 



the dynamics. The point form of relativistic quantum mechanics has important simplifying 



features that are useful in modeling electron scattering. The purpose of this paper is to 
investigate the hadronic current operator in Dirac's point form of dynamics. 

In section II we discuss some of the features of point-form dynamics and construct a 
mass operator for the deuteron. Section III deals with current operators, their relation to 
observables, and the point-form spectator approximation. Then in section IV the numerical 
results are discussed and compared with other methods. Section V presents our conclusions. 

II. POINT-FORM RELATIVISTIC QUANTUM MECHANICS 

Unlike nonrelativistic quantum mechanics, where all the interactions can be put in the 
Hamiltonian operator, for relativistic quantum mechanics it is necessary that at least three 
generators contain interactions. This can already be seen by examining the commutator of 
the Lorentz boost generators with the momentum generators. Such a commutator produces 
the Hamiltonian; if the Hamiltonian contains interactions, then some combination of boost 
and momentum generators must also contain interactions. In the instant form additional 
interactions are put in the boost generators, leaving the momentum generators free of inter- 
actions, while in the point form the additional interactions are in the momentum generators, 
with the boost generators free of interactions. The front form puts interactions in a mixture 
of Lorentz and momentum generators. 

Even though all forms of dynamics are scattering equivalent, each has certain advantages 
that are useful for specific applications. The goal of this paper is to analyze elastic deuteron 
form factors using the point form. The point form has a number of features that set it aside 
from the other forms. First, all of the interactions are in the Hamiltonian and momentum 
generators, that is, the four-momentum operator. Since there are no interactions in the 
boost or angular momentum generators, the Lorentz generators are all kinematic and the 
theory is manifestly Lorentz covariant. It is convenient to write the Poincare commutation 
relations not in terms of the ten generators, but rather in terms of the four-momentum 
operators that contain the interactions, and global kinematic Lorentz transformations: 
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[P M> £,] = 0; (5) 
U A P,U^ = (A-%P„; (6) 

where U\ = U\(A, 0) is a unitary operator representing the Lorentz transformation A. These 
rewritten Poincare relations will be called the point-form equations, and are the fundamental 
equations that have to be satisfied for the system of interest. The mass operator is given by 



M = y P ■ P and must have a spectrum that is bounded from below. 

Since the interactions are all in the four-momentum operators, which are the generators 
of space-time translations, the nonrelativistic Schrodinger equation can be generalized to a 
Lorentz covariant relativistic Schrodinger equation, namely 

idV x /dx» = P^ x , (7) 

where x = x^ is the four- vector space-time point. If the four- momentum operator does not 
depend explicitly on space-time, this equation becomes the eigenvalue equation 

Pfl = P^- (8) 

Finally, as will be shown in the following paragraphs, it is possible to define states with 
the property that angular momentum can be coupled in exactly the same way as is done in 
nonrelativistic quantum mechanics. 

The simplest example of a system satisfying the point form equations is a one-particle 
system with mass m and spin j. If \p, a) is an eigenstate of four-momentum p (with p-p = m 2 ) 
and spin projection er, then 

Pulp,*) =Pi*\p,<r) (9) 

U A \p,a)=^\AprfI%„(R w )M, (10) 

V 

with Ry/ a Wigner rotation defined by R\y = B~ 1 (Av)AB(v), and B(v) a canonical spin 
(rotationless) boost (see reference P?f| ) with argument v = p/m. D ] ala .(R w ) is a Wigner D 
function, and the eigenstates are normalized to 
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(p',a'\p, ( r) = 5- i (p'-p)5 a 



(11) 



relativity requiring the \Jv' q /vq factor. 

States of many noninteracting particles are tensor products of one-particle states; how- 
ever a problem arises when such many-particle states are Lorentz transformed. As can be 
seen from Eq. (jT^) each state is Lorentz transformed by its own Wigner rotation, which in 
general are different. This means that these multiparticle states cannot be directly coupled 
together as is the case nonrelativistically. Such a problem is resolved by coupling the single- 
particle states in the overall rest frame and boosting. It is convenient to label the state by 
the system's four-velocity v: 



ki,//j) := U B ( v )(\ki, /ij ...\k n , /j, n )) 

= Yl \Pl><Tl)-\Pn><7n) H 



0i)o 



(12) 



where Pi = B(v)ki, ki = 0, and Ryy. = B 1 (p i /m)B(v)B(k i /m). Under Lorentz transfor- 
mations, using the definition, Eq. (|12|)> such velocity states transform as 



U A \v,k h [Xi) = | Aw, iZwki, a4) II 



(13) 



where the Wigner rotation Ry/ = B~ 1 (Av)AB(v) is the same in all the arguments of the D 
functions and all the internal momentum vectors ki. That means all the spins as well as the 
orbital angular momenta can be coupled together exactly as is done nonrelativistically. This 
property will be used in the following paragraphs for coupling the nucleon spins together 
with the relative orbital angular momentum to get the spin of the deuteron. From the 
relation between external and internal momenta, it follows that the velocity states defined 
in Eq. ([12]) are eigenstates of the noninteracting mass operator Mf ree and free four-velocity 
operator V^: 



Mfreel^ki,/^) = ^ \J mf + k^ \v, k u fa) ; 

i 



(14) 
(15) 
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The Bakamjian-Thomas procedure is implemented in the point form by writing = 
MV^, where now M is the sum of free and interacting mass operators, M = Mf ree + Mint- 
M takes the place of the center of momentum Hamiltonian h = H — ^ in nonrelativistic 
quantum mechanics; note however that even though there is only one operator containing 
interactions, namely the mass operator, that nevertheless there are interactions in all four 
components of the four-momentum operator. 

In order that the four- momentum operator satisfy the point-form equations, Eqs. (§,0), 
the interacting mass operator must satisfy certain conditions. To satisfy Eq. (|5]), the mass 
operator must commute with the four- velocity operator, defined in Eq. (|15|): 

M, V"] = 0. (16) 

This has the consequence that mass and four-velocity can be simultaneously diagonalized. 
Eigenstates of the four-momentum operator can thus be written as the mass times the 
four-velocity. Since the four-velocity is purely kinematic, it can be factored from the wave 
function leaving the covariant Schrodinger equation, Eq. (§D, to become a mass operator 
eigenvalue equation, 

M$ = A$. (17) 

Moreover, even though the four-momentum is conserved in reactions, the total four- 
momentum is not the sum of the four-momenta of the individual particles. Rather what 
is conserved is the overall four-velocity of the individual particles, and the mass is then 
"off-shell" , not unlike the situation with Feynman diagrams. This is to be contrasted with 
the instant form, where the three-momentum of all the individual particles give the total 
three-momentum of the system, while the energy is "off-shell" . 

The mass operator must also satisfy the other point form equation, Eq. (|6]), implying 
the mass operator is a Lorentz scalar. On velocity states this means the kernel of the mass 
operator must be rotationally invariant and independent of v 2 , exactly the condition put on 
nonrelativistic Hamiltonians in order that they be Galilei invariant. 
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For a two-body system such as the deuteron, the relevant Hilbert space is the tensor 
product of proton and neutron Hilbert spaces, H = H p ®H n . In that case the velocity states 
can be written as \v, k, /i p ,/x n ), where k = ki = — k 2 , and fi p and fi n are the eigenvalues of 
the three-components of the canonical spins of the proton and neutron respectively. Because 
with velocity states the angular momenta can all be coupled together, these states can also 
be written as \v, \k\,j,fij,l,s), as in the nonrelativistic case. The mass of the two particle 



state, from Eq. fll4f) , is 2\Jm? + k 2 ; j is the total angular momentum, while I and s are the 
orbital and spin angular momentum respectively. 

It is advantageous to express the interacting mass operator in terms of a mass squared 
operator with matrix elements: 

(v, \k\,j,fij,l,8\M]\v', |k'|, j',/4,/',s') 

= 5{v - i/)<W;M*> h s|IK) 2 ll^ l', s). (18) 

A mass operator with a kernel of the form Eq. will satisfy Eq. (|T6|) and thus the Poincare 
commutation relations, Eqs. (§,0). The kernel of Mf is taken to be 

(k, I, s\\{m\) 2 \\k\ I', s') := Am(k, /, s^Jk', /', s'), (19) 

where v 3 nn is a nucleon-nucleon interaction. The mass is then defined by 



.\/ : \ M-: M := 4(k + m ) + Mf . (20) 
Denoting the eigenvalue of the interacting mass operator by A 2 , the equation 

M 2 $ = (4m 2 + 4k 2 + 4m< n )$ = A 2 $ (21) 



can be rewritten |24| in the form of the nonrelativistic Schrodinger equation, 



'£ + •*•)•= (£-")•■ (22) 



This defines a relativistic model of the two-nucleon system. It can be shown ^2 that this 
model leads to a small correction to the nonrelativistic binding energy and has scattering 
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observables identical to the corresponding nonrelativistic model. Equation ([TJ) shows that 
the solution of Eq. ( f22"|) leads to simultaneous eigenstates of the mass, velocity, spin, and 
z-component of spin. The Poincare transformation properties of the deuteron eigenstates 
are given by 

P^v^nJiHj) = Xv^v^dJ,^) (23) 

and 



j 3 V v o 

where 

(v, |k|, j,fij,l, 8\v',m D ij\n'j) = 5(v - v')5^ N 5 rj ^\ s {\k\). (25) 

^j a (|k|) is the nonrelativistic wavefunction associated with one of the two chosen nonrela- 
tivistic potentials. (There are analogous formulas for the scattering states.) This provides 
the desired point-form dynamics. 



III. CURRENT OPERATORS, FORM FACTORS, AND ELASTIC OBSERVABLES 



The second key element in a theoretical description of electron-scattering is a conserved, 
covariant hadronic current density J M (x). In the point form the dynamical Poincare trans- 
formations are the space-time translations. Translational covariance can be realized by using 
the dynamical translation operators to define J^(x) in terms of J M (0): 

JP(x) := e ip - x J»{Q)e- ip x . (26) 



The density J^(0) is assumed to transform as a four- vector |p~4"| , p8| with respect to the free 
Lorentz transformations. 

We now want to define the model current operator in terms of measured one-body current 
operators. This is done as follows. The deuteron matrix elements of J M (0) are defined in 
terms of their Breit-frame values with Q in the z direction: 
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(Q/2,l,/4|>(0)|-Q/2,l,^>. (27) 

For n = 0,1,2 the current matrix elements are defined in terms of the single-nucleon 
current matrix elements: 

(Q/2 J l,Ai}|>(0)|-Q/2 ) l ) Ai i > = 

(Q/2, 1,^.| (j£(0) ® I n + I P ® J£(0)) | - Q/2, 1,/i,). (28) 
Current conservation requires that 

]T Q M (Q/2, 1, /4|>(0) | - Q/2, 1, ^) = 0, (29) 
which generates a dynamical contribution Jp n (0) to the z component of the current: 
(Q/2,1,^|J^(0)|-Q/2,1,^-) 

= -(Q/2, 1, (J p 3 (0) <g> J n + J p ® J3(0)) | - Q/2, 1, //,-). (30) 

These relations define the components of the Breit-frame matrix elements of J M (0). The 
remaining deuteron matrix elements of J^(x) are fixed by kinematic Lorentz covariance 
and dynamical space-time translational covariance. Although the current matrix element is 
defined in the Breit frame, the expression for the general current matrix element is Lorentz 
covariant, as can be seen in reference |[28|| , Eq. 3.31. 

The computation of the matrix elements is carried out by inserting single-particle inter- 
mediate states in the velocity basis that was used to formulate the dynamical model in the 
previous section. The deuteron wavefunction in the basis Eq. (p5|) has the form 



(v, \k\,j,fj, j ,l,s\v',m D ,j,tJ l ' j ) = S(v - v^S^^Sj i[Si u (k) + 5 l2 u 2 {k)}, (31) 

where uo(k) and u 2 {k) are the nonrelativistic S and D state deuteron wavef unctions. Trans- 
formation coefficients [p2| , [27| are used to express this in terms of single-particle basis ele- 



ments: 
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(Vi,fl 1 ,V2,fi2\v',fl' i ,m D ) 



5 3 [v / -v(v 1 ,v 2 )] 



^[fc-fc( Vl ,v 2 )] 
k 2 



«9(v,k) 



d(vi, v 2 ) 



1/2 



><<K[ 5_1 K) 5 (^) 5 ( fc i)]<^ 5_1 (^)5M5(A: 2 )]^[ki(v 1 ,v 2 )] 
x(-,//i, -,h 2 \1,/j, s )(1,iii, {<y/ «o[A;(vi, v 2 )] + 5j 2 w 2 [M v i> v 2 )]} 



(32) 



These expressions can be used to compute the current matrix element 

(v,m D , l,[Mj\jg A (0)\v',m D , 

(v, m D , 1, fij\vi,m,V2, fJ-2)(v[, v' 2 , /i' 2 \v', m D , 1, ^J) 



X 



Ml 



(33) 



where the nucleon current matrix elements are given by Eq. (E8) . After integrating out the 
delta functions, one is left with a final three-dimensional integral: 

{v\m D X^\Hu)\v,m d ,l, N ) = E E J d 3 k 

pi 1/2 1/2 wl 1/2 1/2 pin plJ'l 

x^;,(^,0'K(|k'|) x F^(0 )( ^(|k|) 
x^tfl^^i^K)]} x ^{^[^,5(10]} 
x£>^{J2Hr[*i,S(t;)]} x iSgjJMfa.flO;)]} 



x^xht^Kpi -pi) 2 ] + ^e^ (p L Pi) ^ 2[(p; -po 2 ]}^!^; 



2m 



iV 



-{1^2}; 



(34) 



where the C's are Clebsch-Gordan coefficients and the conventions for the spinors, gamma 



and sigma matrices are those of Bjorken and Drell [29 



In this form it can be seen that the momentum of the unstruck particle (the spectator) 
is unchanged, while the struck particle's momentum is changed, but the impulse given to 
the struck particle is not the impulse given to the deuteron. For this reason we call this 
the point-form spectator approximation (PFSA). It should not be confused with the use of 
the term spectator approximation in, for example, reference |3J]. The practical advantage of 
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the PFSA is that the steps above can be generalized to any hadronic target. Moreover, the 
current matrix element is generally Lorentz covariant and can be evaluated in any frame. 

Because the interactions in the point form are in the four-momentum, in the PFSA the 
momentum transfer seen by the scattered nucleon is not the same as the momentum transfer 
seen by the nucleus. In Appendix A we show that the relationship between the momentum 
Q transferred to the deuteron and the momentum transferred to the interacting nucleon is 



Ik"*) l = G— ^T— (1 + 4^) ( 35 ) 



4m|r n _Q 2 , ^ ^ 2 
m D Am D 



>Q 2 ^(l + -^-)>Q\ (36) 



That is, the point-form momentum transfer seen by an individual nucleon will be greater in 
magnitude than the total deuteron momentum transfer Q 2 . 

Two important implications follow from the equations (|35|) and (p6|) above. First, the 
PFSA momentum transfer depends on the internal momentum k, which is a variable of 
integration. Thus in the PFSA, form factors depending on {jp[ — pi) 2 must remain inside the 
integral. Second, since \{p[ — Pi) 2 \ > Q 2 , the deuteron form factors will fall off faster in the 
point-form calculations than in forms where \{p[ — Pi) 2 \ = Q 2 ■ 

The input to the PFSA are single-nucleon current operators. The general structure of 
these operators follows from covariance, parity, hermiticity, and time-reversal symmetry. For 
a spin-1/2 target the conditions imply that all matrix elements can be expressed in terms 



of the Dirac form factors, Fi(Q 2 ) and F 2 {Q 2 ). The general expression has the form [29 



(p,v\J»(0)\p', u')=u u ,(p') 



Fi (q 2 )y + f 2 (q 2 )y: 



u„(p), (37) 



2m 

where u u (p) and u u '{p') are Dirac spinors. In this form the one-body matrix elements are 
easily evaluated in any kinematic frame. 

The Sachs electric and magnetic form factors of the nucleons are 

G E (Q 2 ) = v / TT7(Q/2, 1, l -\ J°(0)| - Q/2, i I)- (38) 

Gm(Q 2 ) = \/^(Q/2, \ \\J\n ~ Q/2, i -\)\ (39) 
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where r = Q 2 /4m 2 . (Here the standard frame is the Breit frame, where the nucleon enters 
with momentum —Q/2 and exits with momentum Q/2, both along the z-axis, which is also 
the axis along which the spin projection is measured.) These Dirac and Sachs form factors 
are related by 

Fi(Q 2 ) = ^-\G E {Q 2 ) + tG m (Q% (40) 

1 + T 

F 2 (Q 2 ) = ^-[G M {Q 2 ) - G E {Q 2 )}. (41) 

1 + T 

The input we use to define the model PFSA current are the single-nucleon form factor 
parameterizations of Gari-Krumpelmann |30| and Mergell-Meissner-Drechsel [ |3T|j . 



The experimental observables for the deuteron and nucleon are well known. The elastic 
observables A(Q 2 ) and B(Q 2 ) are extracted from the Rosenbluth formula for the cross- 
section of unpolarized scattering in the lab frame, 

where a is the fine-structure constant, 9 the scattering angle, and E and E' the initial and 



final energies. For the nucleons, it can be shown |32] that 



Am = gMglilgMgl (43) 

1 + T 

B(Q 2 ) = 2tG 2 m (Q 2 )- (44) 

For spin- 1/2 particles, measurements of A(Q 2 ) and B(Q 2 ) suffice to determine G E and Gm- 
Various models of the nucleon form factors [p0l^lj^3| -|35| have been constructed. 

The deuteron has three independent form factors. A common classification is to denote 
them as the charge monopole G E , magnetic dipole Gm, and electric quadrupole Gq form 
factors. As current matrix elements, these are defined in the Breit frame: 

G i; = i(Q/2,l,0|J (0)|-Q/2,l,0) 

+ |(Q/2,1,1|J°(0)|-Q/2,1,1); (45) 

G M = y^(Q/2, 1, 11^(0)1 - Q/2, 1, 0); (46) 
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G Q = ^[(Q/2,1,0|J°(0)|-Q/2,1,0> 

- (Q/2,1,1|J°(0)|-Q/2,1,1)]; (47) 



where r\ = Q 2 /4m 2 D . These form factors have the static limits 



G E (0) = e; (48) 

lim G M (Q 2 )=e^fi D ; (49) 
Q 2 ->o thin 

YimG Q {Q 2 )=em 2 D Q D] (50) 

where e is the charge, /in the magnetic dipole moment, and Qe> the electric quadrupole 
moment of the deuteron. 

The Rosenbluth formula alone cannot determine all three of the deuteron's form factors. 
The other independent observable normally measured is the deuteron tensor polarization 
T 20 , defined as 

T 20 := V2 d -^^, (51) 

where da^ refers to the differential cross-section with helicity /i. Conventionally it is dis- 
played at a 70° angle in the lab frame. The deuteron elastic observables are: 

A(Q 2 ) =G 2 E + 8 - V 2 G 2 Q + 2 - V G 2 M , (52) 
B(Q 2 ) = ^v(l + v)G 2 M ; (53) 

T 20 (Q ) - -^V Am + Bmtan2{d/2y (54) 

where / = 1/2 + (1 + n) tan 2 (^/2). 

IV. NUMERICAL RESULTS AND COMPARISONS 

One purpose of this work was to test the point-form spectator approximation on the 
simplest nucleus, the deuteron, where realistic interactions and nucleon form factors are 
available. Comparisons with other models are then an indication of the relative size of the 
required two-body currents. 
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For the nucleon-nucleon interaction, we have used the Argonne v±s [[36] and Reid '93 [37 



potentials to construct a mass operator M. The S (1=0) and D (1=2) state wavefunctions 
are plotted in Figures [1] and 0. The only significant differences between the wavefunctions 
these potentials produce in configuration space occur below 0.4 fm for the S wave and below 
1.0 fm for the D wave. In momentum space the wavefunctions do not exhibit significant 
differences up to 5 fm -1 , about 1 GeV, above which they do differ noticeably. The effects 
on the choice of interaction may therefore be expected to be relevant at higher momentum 
transfers, but as will be seen, these high-momentum differences in the wavefunction make 
only slight differences in the calculations. 

The PFSA currents are constructed using the Gari-Krumpelmann [3^] and Mergell- 



Meissner-Drechsel ]3T| parameterizations of the nucleon form factors. At the range of mo- 



mentum transfer under consideration, the parameterizations give very similar results for the 
proton form factors and the neutron magnetic form factor. The neutron electric form factor, 
however, varies significantly between the two. 

The deuteron form factor Ge has been calculated using both form factor parameteriza- 
tions and both nucleon-nucleon potentials. The absolute values of the results are displayed in 
Figure |3|. The results are independent of the nucleon-nucleon potential used, except for small 
variations at high momenta. The primary differences in Ge and Gq are due to the nucleon 
form factor parameterizations. (Figure |4] compares the G-K and MMD neutron form factors. 
Note that the major difference is in the parameterization of the neutron electric form factor; 
the neutron magnetic form factor parameterizations, and the proton parameterizations as 
well, are very similar.) For Ge, both the Gari-Krumpelmann and the Mergell-Meissner- 
Drechsel form factors agree at low momentum transfers and have zeroes near 0.8 GeV 2 . The 
G-K form factors predict a second zero near 5.5 GeV 2 while MMD predicts a second zero 
between 6 and 7 GeV 2 . Because the two form factors are almost identical except for the 
parameterization of the neutron electric form factor, this would suggest that the neutron 
form factor is the dominant cause of the differences in the calculations of Ge- 

Figure || illustrates the dependence of the magnitude of the form factor Gm on the 
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potential and on the nucleon form factors used. Both parameterizations predict the same 
behavior up to the first zero, this time at 1.6 GeV 2 , and within the range studied, fall off 
with almost identical behavior. Comparison with Figure [| would suggest that the neutron 
electric form factor has little effect on the calculation of Gm- A further comparison to 
experimental data can be made by examining the static limit of Gm- Equation (EUp relates 



this limit to the deuteron magnetic moment, and Table H displays the results. In the static 
limit, the parameterization of the nucleon form factors does not affect the results, while 
the choice of nucleon-nucleon interaction does. This is expected, as the form factors must 
approach precise limits as Q 2 — > 0, while the momentum-space wavefunctions have no such 
constraints. 

This procedure is repeated for Gq in Figure ||. As was the case with Ge, there is little dif- 
ference due to the potential used, but a noticeable difference between the predictions of the 
G-K and MMD parameterizations. The G-K form factors show a zero between 4.5 and 5.0 
GeV 2 , while the MMD form factors produce a zero approximately one GeV 2 higher. Further, 
the magnitude of the G-K results is greater than that of MMD almost everywhere through- 
out. The different neutron electric form factor parameterizations is the primary cause of the 
differences in the results. The deuteron electric quadrupole moment (Eq. |5U| ), displayed in 
Table 1, differs from the experimental result, the calculated values approaching about 90% 
of the experimental value, as opposed to 99% for the magnetic moment calculations. This 
is consistent with other models pT| , p!7| , |34| , |36| , |38| -|40[ . 



To summarize, these point-form calculations imply that the deuteron form factors are 
essentially independent of the potential (Argonne Vi 8 or Reid '93) used, but depend more 
significantly on the parameterizations of the form factors, and in particular on the neutron 
electric form factor, as this is the only substantial difference between the G-K and MMD 
parameterizations. The static moments are similar to predictions in other realistic models 
with the electric quadrupole moment differing with experiment by about 10%. 

Figure [7] displays the results for A(Q 2 ) up to 2 GeV 2 for both potentials and both form 
factor parameterizations; Figure [8] extends the calculations to 8 GeV 2 . The data come from 
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Refs. fil|- fP]l . Differences among the various calculations begin to appear at intermediate 
momentum transfers. For A(Q 2 ), with Q 2 between 0.5 and 3 GeV 2 , the PFSA combined with 
the Gari-Kriimpelmann form factors fit the data fairly closely, while the Mergell-Meissner- 
Drechsel form factors produce results that fall short by an order of magnitude. This pattern 
occurs in other impulse and spectator calculations as well. In the front-form calculations of 
Chung et al. |1| , the fit for an earlier G-K parameterization (using the Argonne vu poten- 



tial) is even better, while the Hohler (on which the newer MMD form factors were based) 
calculations again fall an order of magnitude short. In Lev et al. ]ll| the two curves are 



closer, though the G-K remains higher and fits the data out to 2.0 GeV 2 . In the nonrel- 
ativistic calculations of Carlson and Schiavilla [39], which cover the range 0-2.4 GeV 2 and 



use only the Hohler form factors, impulse approximations using various potentials (including 
the Argonne Dig) all fall nearly an order of magnitude short in the intermediate range. In 
the work of Van Orden et al. [[?[, which contains an impulse approximation that fits most of 
the data for all three form factors quite closely, the variation from the A(Q 2 ) data starts at 
2.0 GeV 2 and is an order of magnitude short at high (8 GeV 2 ) momentum transfers. 

The elastic observable B(Q 2 ), related directly to the magnetic form factor Gm, is dis- 
played in Figure up to 8 GeV 2 using both potentials and both sets of parameterizations. 



The data come from Refs. [41J45-47H. The PFSA calculation of B(Q 2 ) in the intermediate 
region 0.5-3 GeV 2 , a region which contains all presently available data, fits that data poorly, 
though the differences between the two form factor parameterizations are less marked. Both 
parameterizations exhibit a zero at 1.5±0.1 GeV 2 , causing a wide discrepancy with experi- 
ment; data suggest a zero nearer 1.9 GeV 2 . Again Chung et al. and Van Orden et al. fit the 
B(Q 2 ) data quite well. In contrast, Chung et al, using the Paris and the Bonn wavefunctions 
instead of the Argonne t>i4, produce results quite similar to the point form's. Carlson and 
Schiavilla obtain a zero at 2.2 GeV 2 only with the Nijmegen potential; their other potentials 
reproduce the zero at 1.6 GeV 2 . Lev et al. produce zeroes between 1-2 GeV 2 in the 0-4 
GeV 2 region using various potentials and parameterizations. The position of the zero in all 
these forms seems to be the most salient feature of calculations of B(Q 2 ), affecting as it does 



the deviation from data in the 1-2 GeV 2 range. 

The tensor polarization, T 20 (Q 2 ), is displayed in Figure [0], with data from Refs. P^-|53|. 
As with B(Q 2 ), measurements of T 2 o(Q 2 ) have only investigated the low (0-0.5 GeV) and 
intermediate (0.5-1.5 GeV) ranges of momentum transfer. In both ranges, the MMD and 
G-K parameterizations produce identical results in the PFSA. In the intermediate range, the 
results fall slightly below the data. The impulse approximations of Carlson and Schiavilla 
as well as Lev et al. do this too, while Van Orden et al. and Chung et al. produce curves 
that fit modern data quite closely. 

Finally, it is instructive to compare the results obtained in point-form dynamics to the 
results of nonrelativistic impulse calculations to get some idea of the nature of the relativistic 
effects; and to compare the results in the point form to the same relativistic calculations 
done assuming \(p[ — Pi) 2 \ = Q 2 (that is, pulling the nucleon form factors outside of the 
integral) to examine how the point-form momentum transfer affects the results. 

The results (using the Argonne v is potential and the Gari-Krumpelmann form factors) 
for A(Q 2 ) are displayed in Figure |TT[ For low momentum transfers, all three agree. In 
the range of 1.0-5.0 GeV 2 , the nonrelativistic and PFSA calculations decrease similarly, 
the nonrelativistic curve consistently higher. The constant-Q 2 calculation in this region 
gradually rises from the PFSA to the nonrelativistic curve. Above 5.0 GeV 2 , all three curves 
systematically decrease, nonrelativistic above constant-Q 2 above PFSA. While a relativistic 
treatment is needed as a matter of principle at high momentum transfers, it is clear that in 
these calculations the effects of combining the PFSA with point-form quantum models has 
a tendency to reduce the structure function A(Q 2 ) at high Q 2 . One clear cause of this is 
that 

M-Pi) 2 \>Q 2 ; (55) 

the magnitude of the point-form momentum transfer is greater than the magnitude of the 
nonrelativistic momentum transfer. Because the form factors depend on the magnitude of 
the momentum transfer, they therefore drop off more quickly in the point form. This reduces 
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the point-form results in comparison to the constant-Q 2 calculations, as Figure |TT| shows: 
as the momentum transfer increases, the two calculations diverge further from each other. 

The point-form and nonrelativistic results for B(Q 2 ) and T2o{Q 2 ) do not exhibit as 
dramatic differences as they did for A{Q 2 ). Again in the graphs of B{Q 2 ) and T 2 o(<5 2 ), 
Figures [T2| and [TB|, one sees the similar but increasingly divergent results obtained from the 
constant-Q 2 and the PFSA methods. 



V. CONCLUSION 

This work has used the point form of relativistic dynamics to calculate elastic deuteron 
form factors. The point form stands somewhat between the covariant approaches and direct- 
interaction approaches mentioned in the introduction in that it is on the one hand manifestly 
covariant (because the Lorentz generators are kinematic) and it is the mass that is "off-shell" 
(rather than the energy as is the case with the instant form.) On the other hand the point 
form is one of the forms of dynamics listed by Dirac, in which all of the interactions are in 
the four-momentum generators. Moreover there is a natural way in which one-body currents 
can be introduced in the point form (called the point-form spectator approximation) that 
satisfies the correct Poincare and charge conservation properties. 

We have shown that the PFSA produces results consistent with other impulse and spec- 
tator approximations. Within the range Schiavilla and Riska examined, for example, their 
impulse approximation and the PFSA (using G-K form factors) predict nearly identical 
results for A{Q 2 ) and T 2 o(Q 2 ); and though the zero they predict for B{Q 2 ) falls near 2.0 
GeV 2 rather than the PFSA's 1.6, the fall-off from the data begins near 0.5 GeV 2 in both. 
The calculations of Kobushkin and Syamtomov (before using their approach of reduced 
transition amplitudes) and the PFSA (G-K) results for A(Q 2 ) and B(Q 2 ) nearly duplicate 
each other, as do the results for T2o(Q 2 ) except in the high-momentum range, where no data 
are currently available. And although the calculations of Chung et al. with earlier G-K form 
factors and the Argonne v u potential fall quite close to the data for all three observables, 
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their calculations using Hohler [35] form factors and potentials show the same salient points: 
the deviation from data beginning in A(Q 2 ) near 1 GeV 2 and in B(Q 2 ) near 0.5 GeV 2 ; the 
location of the first zero in B(Q 2 ) between 1.5-2.0 GeV 2 ; and the first minimum in T2o(Q 2 ) 
around 0.8 GeV 2 . The work of Van Orden et al. gives results similar to Chung et al, except 
that the high-Q 2 behavior of A(Q 2 ) and B(Q 2 ) is nearly level at 10 -8 , while Chung et al. and 
the PFSA show gradual decreases at 10 -11 for A(Q 2 ) and at 10~ 8 for B(Q 2 ). In contrast, 
Lev et al. calculate that the Hohler form factors produce results for A(Q 2 ) that lie close to 
the data, while the G-K form factors fit A(Q 2 ) up to 2 GeV 2 but produce results increasingly 
higher than the data thereafter. Their results for B(Q 2 ) and T2o(Q 2 ) are similar to those of 
Van Orden et al. and Chung et al. 

This work also addressed the sensitivity of PFSA results to different nucleon-nucleon 
interactions and different parameterizations of the nucleon form factors. In almost every 
instance it was found that the two nucleon-nucleon potentials produced only slight, if any, 
differences in the form factors and elastic observables. This may not be surprising considering 
that the Argonne t> 18 and Reid '93 nucleon-nucleon interactions produce nearly identical 
momentum-space wavefunctions on the momentum scale of interest. 

Much more pronounced were the differences between the Gari-Krumpelmann and the 
Mergell-Meissner-Drechsel parameterizations of the nucleon form factors. As the momentum 
transfers become higher, the two often predict significantly different results. These are most 
notable in the deuteron form factors Ge and Gq, which are sensitive to the neutron electric 
form factor; the G-K parameterization, whose neutron electric form factor falls off markedly 
more rapidly than the MMD, produces deuteron form factor zeroes in the intermediate range 
that occur at higher momentum transfers in the MMD results. That this phenomenon is due 
to the neutron electric form factor is supported by the similarity of results for Gm, where 
the nucleon magnetic form factors dominate the calculations. 

Aside from differences due to varying potentials and nucleon form factors, the conse- 
quences of the point form's nontrivial momentum transfer have also been examined. We 
have shown that the point-form momentum transferred to a nucleon is greater than the Q 2 
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transferred to the deuteron, and that its deviation increases with increasing Q 2 . This results 
in a lowering of the deuteron form factors and elastic scattering observables compared to 
nonrelativistic calculations. The greater magnitude of the point-form momentum transfer 
causes a quicker fall-off of the nucleon form factors, and some deviation from nonrelativistic 
calculations was attributed to this. 

Additional dynamically consistent two-body currents must be added to the PFSA in 
order to bring the calculations into agreement with data at intermediate to high momentum 
transfers. Such currents have not been considered in this preliminary work; but other 
calculations that include dynamical two-body currents (see references |]8| JlllJT5|JT7|] ) suggest 
that the addition of dynamical currents is capable of reconciling the differences between 
various impulse or spectator approximations with data. 

In the point form it is quite easy to interpret Feynman diagrams for nucleon-nucleon 
scattering with the production of a photon as a current matrix element satisfying the re- 
quirements given in Section III. However, while the addition of such current matrix elements 
may produce better agreement with data, it does not provide a systematic procedure for 
constructing two-body currents. What is needed is a procedure for constructing conserved 
currents from one-body currents and the dynamical mass operator. Models based on such 
a procedure are being developed. 



APPENDIX A: POINT-FORM MOMENTUM TRANSFER 



The momentum transfer (p^ — pi) seen by nucleon i can be computed following Ref. [28 
Suppose that the momentum transfer is along the z-axis. 



B{v in ) 



cosh A/2 sinhA/2 ^ 

10 

1 

sinhA/2 cosh A/2 



(Al) 
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B(v out ) = 



(A2) 



cosh A/2 -sinh A/2 ^ 
10 
1 

^ -sinh A/2 cosh A/2 

are the boosts that take the deuteron from the center of momentum frame to the Breit frame 
(where P' tot = — P to t-) For elastic scattering, 



/ 



sinh A/2 = 
The initial energies and momenta are then: 



Q 2 



Am 2 D 



(A3) 



E x = u cosh A/2 + k z sinh A/2; 
Piz = k z cosh A/2 + uo sinh A/2; 
E 2 = u cosh A/2 - k z sinh A/2; 
P2z = — k z cosh A/2 + uo sinh A/2; 



(A4) 



where u and k are center of momentum variables. 

In this notation, k z refers to the relative z-axis momentum of particle one. This conven- 
tion gives rise to the following relations: 



u' = uj cosh A =)= k z sinh A; 
k' z = k z cosh A =F u sinh A; 



(A5) 
(A6) 



where the minus signs are used when particle one is struck, the plus signs when particle 
two is struck. Suppose for illustration that particle one is struck. The final energies and 
momenta will then be 

E[ = cosh 3 A/2 - k z sinh 3 A/2; 
Piz = kz cosh3A/2 — u sinh3A/2; 

E' 2 = E 2 ; p' 2z =p 2z - (A7) 
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Now some hyperbolic trigonometry reveals that 



( p [ - Pl f = 4(£; 2 - cu 2 ) sinh 2 A 



(A8) 



Since 



sinhA = 2W^U/l+ Q2 



4m 2 D y 4m 2 D ' 



(A9) 



and 



k; - oj = k 2 z - ml - k 2 = -(m 2 N + k : 



(A10) 



the resulting Equation 36 is established. 
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TABLES 



Moment 


Units 


Experimental 


Argonne v±s 


Argonne v±s 


Reid '93 


Reid '93 








G-K 


MMD 


G-K 


MMD 


fJ-D 




0.85741 


0.8613 


0.8623 


0.8615 


0.8625 


Qd 


e/GeV 2 


7.3422 


6.6 


6.6 


6.6 


6.6 



TABLE I. The magnetic dipole and electric quadrupole moments, computed as Q 2 — > 0. 
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FIGURES 

FIG. 1. The bound state S (1=0) deuteron wavefunction, using the Argonne (solid) and Reid 
(dashed) potentials. 

FIG. 2. The bound state D (1=2) deuteron wavefunction, using the Argonne (solid) and Reid 
(dashed) potentials. 

FIG. 3. G E (Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. 

FIG. 4. The neutron electric and magnetic form factors of G-K (dotted and dash-dot), and of 
MMD (dashed and solid). 

FIG. 5. Gm{Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. 

FIG. 6. Gq(Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. 

FIG. 7. A(Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. The data come 
from Refs. [41] (squares), [42] (crossed circles), [43] (triangles), [44] (open circles) and [56] (curved 
squares) . 

FIG. 8. A{Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. The legend for 
the data is the same as in Figure 7. 

FIG. 9. B{Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. The data come 
from Refs. [41] (squares), [45] (triangles), [46] (crossed circle), and [47] (open circles). 
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FIG. 10. T 20 (Q 2 ) for the Argonne potential with G-K (solid) and MMD (dashed) form factors, 
and for the Reid potential with G-K (dash-dot) and MMD (dotted) form factors. The data are 
compiled from Refs. [48-55]. 

FIG. 11. Point-form (solid), nonrelativistic (dashed), and constant-Q 2 (dash-dot) results for 
A{Q 2 ) using the Argonne potential and the G-K form factors. 

FIG. 12. Point-form (solid), nonrelativistic (dashed), and constant-Q 2 (dash-dot) results for 
B(Q 2 ) using the Argonne potential and the G-K form factors. 

FIG. 13. Point-form (solid), nonrelativistic (dashed), and constant-Q 2 (dash-dot) results for 
^2o(Q 2 ) using the Argonne potential and the G-K form factors. 
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